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Growth Dynamics of Photoinduced Domains in Two-Dimensional Charge-Ordered 
Conductors Depending on Stabilization Mechanisms 
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Photoinduced melting of horizontal-stripe charge orders in quasi-two-dimensional organic 
conductors 6>-(BEDT-TTF) 2 RbZn(SCN) 4 [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene] 
and a- (BEDT-TTF) 2 13 is investigated theoretically. By numerically solving the time-dependent 
Schrodinger equation, we study the photoinduced dynamics in extended Peierls-Hubbard mod- 
els on anisotropic triangular lattices within the Hartree-Fock approximation. The melting of the 
charge order needs more energy for 6>-(BEDT-TTF) 2 RbZn(SCN) 4 than for a-(BEDT-TTF) 2 I 3 , 
which is a consequence of the larger stabilization energy in #-(BEDT-TTF)2RbZn(SCN)4. Af- 
ter local photoexcitation in the charge ordered states, the growth of a photoinduced domain 
shows anisotropy. In #-(BEDT-TTF)2RbZn(SCN)4, the domain hardly expands to the direc- 
tion perpendicular to the horizontal-stripes. This is because all the molecules on the hole-rich 
stripe are rotated in one direction and those on the hole-poor stripe in the other direction. They 
modulate horizontally connected transfer integrals homogeneously, stabilizing the charge order 
stripe by stripe. In a- (BEDT-TTF) 2I3, lattice distortions locally stabilize the charge order so 
that it is easily weakened by local photoexcitation. The photoinduced domain indeed expands 
in the plane. These results are consistent with recent observation by femtosecond reflection 
spectroscopy. 

KEYWORDS: photoinduced phase transition, charge order, metal-insulator transition 



1. Introduction 

Low-dimensional organic conductors are known to of- 
fer a playground for studying strongly correlated elec- 
tron systems in which various interesting ground states 
emerge. 1 ' 2 The typical examples are Mott and charge- 
ordered insulating states, where electron-electron as well 
as electron-phonon (e-ph) interactions have important 
roles. In such systems, photoirradiation may induce a 
phase transition, which has been extensively studied re- 
cently 3 6 since it will lead to novel transient phenomena 
and possible control of various functions of materials. 

The quasi-two-dimensional organic conductors 6- 
(BEDT-TTF) 2 RbZn(SCN) 4 (abbreviated as (9-RbZn 
hereafter) and a-(BEDT-TTF) 2 I 3 (as a-l 3 ) are typical 
compounds that exhibit a charge order (CO). 7 10 They 
consist of stacking layers of monovalent anions and donor 
BEDT-TTF molecules whose 7r-band is 3/4-filled. The 
photoinduced melting of the COs in these systems has 
recently been observed by using femtosecond reflection 
spectroscopy. 6, 11 It shows a marked difference between 
their photoinduced dynamics. For a-ls^ a semimacro- 
scopic metallic domain is generated, whereas the CO only 
locally melts for #-RbZn. In particular, the dynamics in 
a-ls shows critical slowing down and strong dependence 
on the excitation intensity and temperature. 6 

This clear difference is considered to originate from dif- 
ferent roles of lattice distortions in stabilizing the COs. 
In fact, the CO transition in #-RbZn is a first-order 
metal-insulator transition with large structural distor- 
tion at T c = 200K. 12,13 On the other hand, a-ls under- 



*yasuhiro@ims. ac.jp 
tkxy@ims.ac.jp 



goes small lattice distortion at the CO transition with 
T c = 135K. 14,15 In both salts, a horizontal-stripe CO 
is formed, as confirmed by several experiments such as 
X-ray scattering 13, 15 and Raman spectroscopy. 16, 17 

Theoretically, CO phenomena have been investigated 
by using extended Hubbard models including on-site (U) 
and intersite (Vij) Coulomb interactions. 18 35 The stabil- 
ity of various CO patterns in some members of (BEDT- 
TTF) 2 X has been discussed first within the Hartree ap- 
proximation by considering realistic band structures of 
(BEDT-TTF) 2 X. 1,18 The horizontal-stripe CO is shown 
to be stabilized in #-RbZn and also in a-l3 18 if the trans- 
fer integrals are based on the low-temperature crystal 
structures with lattice distortions implicitly included. 
The effects of e-ph interactions on the horizontal COs 
have been investigated by the Hartree-Fock approxima- 
tion and exact diagonalization. 30 34 Here the transfer in- 
tegrals are based on the high-temperature crystal struc- 
tures if the lattice is undistorted. The results show that 
charge frustration due to the Coulomb interactions on a 
triangular lattice and lattice distortion relieving the frus- 
tration are essential for #-RbZn, while the low-symmetry 
configuration of transfer integrals is important for a-ls. 
These theoretical studies consistently interpret the ex- 
perimental findings including the large discontinuity at 
the CO transition in #-RbZn and the difference between 
the behaviors of the spin degrees of freedom in the two 
salts below T c . 12,36 Therefore, it is of great importance 
to investigate the photoinduced melting of the COs and 
to compare their dynamics. 

In this paper, we investigate the photoinduced melt- 
ing dynamics in #-RbZn and in a-13 using the time- 
dependent Schrodinger equation for extended Hubbard 
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models with Peierls-type e-ph couplings within the 
Hartree-Fock approximation. It is found that the CO 
in #-RbZn is more stable against photoexcitation than 
that in a-ls, which is consistent with the experimental 
results. 6 Their dynamics show different behaviors owing 
to different roles of e-ph couplings for the COs and dif- 
ferent underlying crystal structures. In particular, a pho- 
toinduced domain with suppressed CO grows anisotrop- 
ically in #-RbZn, whereas it grows isotropically in a-13. 
This suggests that a macroscopic domain is more eas- 
ily created in 0-I3 than in #-RbZn. In §2, the extended 
Peierls-Hubbard model is defined. The numerical method 
for solving the time-dependent Schrodinger equation is 
also given. After presenting linear absorption spectra in 
§3, we show photoinduced dynamics during and after a 
spatially uniform, oscillating electric field is introduced 
in §4. The growth of a photoinduced domain after local 
melting of the CO is discussed in §5. Section 6 is devoted 
to a summary. 

2. Extended Peierls-Hubbard Model on Trian- 
gular Lattice 

In order to describe the horizontal CO with lattice 
distortion in #-RbZn and a-ls, we consider the following 
extended Peierls-Hubbard model 30 34 

H = H e \ + H\ at , (1) 

with 

+ u E n *r n a + E V hj n i n j ' 

(ij) (ij) 

where (ij) represents the summation over pairs of neigh- 
boring sites, cl a (ci a ) denotes the creation (annihilation) 
operator for an electron with spin a at the ith site, 
nia = c \a c i(j, and Tti = ritf + 71^. The electron density is 
fixed at 3/4 filling. The on-site repulsion is denoted by U. 
For the inter site Coulomb interactions Vij, we consider 
nearest-neighbor interactions V c in the vertical direction 
and V p in the diagonal direction, as shown in Fig. 1(a). 
The e-ph coupling constant, lattice displacement, elas- 
tic constant and bare phonon frequency are denoted by 
otij, Uij, Kij and ojij, respectively. We introduce new 
variables by jjij = ctijUij and Sij = afj/Kij. The no- 
tations of the transfer integrals Uj, yij and Sij are the 
same as those in refs. 32 and 33 for both compounds. 

The structures of #-RbZn and a-ls in the high- and 
low-temperature phases are schematically shown in Fig. 
1. In the metallic phase of #-RbZn, there are two kinds 
of transfer integrals t p and t c for the diagonal and verti- 
cal bonds, respectively, while the CO transition doubles 
the unit cell in the c-direction and six transfer integrals 
appear at low temperatures. On the other hand, the unit 
cell of a-ls contains four molecules in both phases. Ac- 
cording to the X-ray structural analysis, 15 sites A and 
A' are equivalent in the metallic phase owing to inver- 



sion symmetry, while the symmetry breaks below the CO 
transition temperature. 

In the previous studies, 30 ' 32 ' 33 we consider three types 
of e-ph couplings s c , s p i, and for #-RbZn, and two 
types of e-ph couplings Sbi and for a-ls in order to 
investigate the roles of all the lattice distortions on the 
stability of the COs. For #-RbZn, the high-temperature 
structure has high symmetry, where several patterns of 
COs are nearly degenerate because of the charge frustra- 
tion. The lattice distortions relieve the frustration. We 
found that s^, which comes from molecular rotation, 30 ' 31 
is the most important for realizing the horizontal CO. As 
shown in Fig. 1(c), this distortion homogeneously modu- 
lates the horizontally connected transfer integrals, 30 ' 31 
t P 2 = tp — Vcf) and t p 4 = t p + 2/0, by which the CO 
is globally stabilized. For a-ls, on the other hand, lat- 
tice effects are relatively small and the low-symmetry 
configuration of transfer integrals stabilizes the horizon- 
tal CO, 32 ' 34 as explained below. Charge disproportion- 
ation between sites B and C exists even in the metallic 
phase 38 41 because the large and nonequivalent transfer 
integrals t^i and £52 form a zigzag chain, 15 as seen in Fig. 
1(b). Since is larger than t^i in the chain, site B (C) 
becomes hole-rich (hole-poor) so as to gain the kinetic 
energy. The appearance of the horizontal CO breaks the 
equivalence of sites A and A' . The e-ph couplings s^i 
and 552 stabilize this CO. However, the configuration of 
large transfer integrals is not along the hole-rich stripes 
In contrast to the case of #-RbZn. In other words, the 
'CO is locally stabilized in a-ls: the distortion causing 
tB2' > ts2 (tBv > tsi) in Fig. 1(d) locally strengthens 
(^ie bond B2' (Bl 7 ), which is regarded as the formation 
of a local singlet on the bond B2' . As a consequence, site 
A (A') becomes hole-rich (hole-poor). With these facts 

mind, in the present paper we consider only and 
as the e-ph coupling in #-RbZn and a-ls, for simplicity. 
Because only one phonon mode is taken into account, we 
set the bare phonon frequency at Uij = c^ph- 

The photoexcitation is introduced by the Peierls phase 
factors of the transfer integrals in eq. (2). Here Sij de- 
notes the position vector from the ith site to the jth site. 
The vector potential A(t) is given by, 

A(t) = -c [ dt'E(t') , (4) 
Jo 

E(t) = E ext 0(t)0(T [rr - t) sinu ext t , (5) 

where 22 e xt and Cc; ex t is the amplitude and frequency, re- 
spectively, of the time-dependent electric field E(t). In 
eq. (5), we define the pulse width T- 1TT = 27r7V ext /cj ext 
with iVext being an integer. 0{t) is the Heaviside step 
function, 0(t) = 1 for t > and 9(t) = for t < 0. We 
use e = 1, h = 1, and the intermolecular distance along 
the c-axis of #-RbZn and that along the a-axis of a-ls 
as the unit of length. Furthermore, the intermolecular 
distance along the a-axis of #-RbZn and that along the 
6-axis of a-ls are set at twice the unit length. By these 
definitions, parameters are given in units of eV unless 
otherwise noted. 

The time evolution of the system is calculated as fol- 
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t a j = -3.5(xl(T 2 eV) 
t.2 = -4.6 
t*s = 1-8 
tbi = 12.7 
t b2 = 14.5 

tbs = 6-2 
t b4 = 2.5 




td = 1.5 
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Fig. 1. (Color online) Schematic representations of the structures 
of (a) <9-(BEDT-TTF) 2 X in the metallic phase, (b) a-h in the 
metallic phase, 37 (c) #-RbZn in the CO phase, and (d) CK-I3 in the 
CO phase. The gray solid lines indicate unit cells. For (b), (c), 
and (d), the transfer integrals estimated by the extended Hiickel 
method 13 ' 15 are also shown. The red or thick solid (dashed) lines 
in (c) and (d) indicate the bonds on which the magnitudes of the 
transfer integrals are increased (decreased) by the distortions 
considered in the present paper. 



lows. For the electronic part, we use the time-dependent 
Schrodinger equation, 

\<i/j^ a (t + At)) =Texp 



dt>H* F (t>) |W(*)>, 



(6) 

where T denotes the time-ordering operator, \i/j u ^(t)) the 
v-th one-particle state with spin a at time t, and H^ F the 
electronic part H e \ in the Hartree-Fock approximation. 

For the lattice part, we solve the classical equation 
of motion by the leapfrog method. 42 In the case of 0- 
RbZn, we assume that for the lih unit cell the transfer 
integrals on the two p4-bonds connected to site B and 
those on the two p2-bonds connected to site C are given 
by t p4 (l) = t p + 2/0(0 and t p2 (l) = t p -y<j,(l), respectively. 
The equation of motion for u<f>(l) is then written as, 

(fu^il) 



dt 2 



1*7®, 



where we write the force for #-RbZn as, 



#<!) = -(- 



*(*)> 



(7) 



(8) 



-du^l) du<f>(l) 

Here = K^juJK is the phonon mass, and is the 

Slater-determinant composed of the one-particle states 
\ipv i(T (i)). In the case of 0-I3, we set t B 2'{l) = h 2 +2/62(0 
and ^2(0 = tb2 — 2/62(0- The equation of motion for 
1^2 (0 thus reads, 



m 6 2- 



d 2 u h2 (l) 
dt 2 



= - 2 mi 



(9) 



where = if 62 /^ph anc ^ * s defined as in eq. (8). 

We note that for both compounds the initial state is the 
horizontal CO state where the lattice displacements do 
not depend on Z. 30-34 When the uniform time-dependent 
electric field is applied, their response is still uniform. 
However, they become nonuniform when the electric field 
is applied locally as in §5. 

In the actual calculations, we first obtain Ff(t) or 
F t a (t) by using the wave function \ifj u ^(t)). Then, or 
2/62 at t + At is calculated, which is accurate to the order 
of (At) 2 . The wave function |Vv,cr(^ + At)) is obtained 
from eq. (6) by, 



^(t+At)) ~exp 



1 



(10) 



whose error is of the order of (At) 3 . In the above equa- 
tion, we use H* F (t + \ At) = [H* F (t) + H* F (t + At)]/2 
where H^ F (t + At) is constructed so that it is accu- 
rate to the first order of At. The exponential operator 
is expanded with time slice At = 0.01 until the norm of 
the wave function becomes unity with sufficient accuracy. 
The obtained \ip u< r(t + At)) has an error of the order of 
(At) 3 . 

In this paper, the phonon frequency c^ p h = 0.01 is used 
for both compounds. We set t p = 0.1 and t c = —0.04 
for #-RbZn. For the transfer integrals in a-13, we use 
the values shown in Fig. 1(b). The system size is fixed 
at 12 x 12 sites and the periodic boundary condition is 
imposed. The interaction strengths and the e-ph coupling 
constants are chosen at U = 0.7, V c = 0.4C7, V p /V c = 0.6, 
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Fig. 2. (Color online) Linear absorption spectra obtained by the 
time-dependent Hartree-Fock approximation for 12xl2-site sys- 
tem with £7=0.7, V c /U=0A, and V p /V c =0.6 for (a) <9-RbZn 
(s0=O.l) and (b) ai-13 (s52=0.07) with different polarizations 
as indicated. The peak-broadening parameter 7 is set at 0.02. 



= 0.1, and S52 = 0.07 throughout the study. These 
parameters give the horizontal CO ground states for both 
compounds. 30 ' 32 

3. Linear Absorption Spectra 

First, we discuss linear absorption spectra, which cor- 
respond to optical conductivity spectra. For this purpose, 
we replace the electric field in eq. (5) by the one with a 
damping factor, 



E(t) = E ext 6(t) exp(-7t) sincjext^ 



(11) 



and calculate the increment in the total energy due to 
E(t) with small | E ext | at sufficiently large t so that 
the total energy converges. We use | i£ ext | =0.002 and 
7=0.02 to obtain results shown in Fig. 2. In both cases 
of #-RbZn and the absorption is larger for the polar- 
ization parallel to the stripes, i.e., along the a- and 6-axes, 
respectively, than for the polarization perpendicular to 
them, i.e., along the c- and a-axes, respectively. This is 
because the transfer integrals on the diagonal bonds are 
larger than those on the vertical bonds. The absorption 
spectra with polarization parallel to the stripes have sim- 
ilar structures with resonance peaks located at around 
uj = 0.4 for both compounds. Therefore, in the follow- 
ing we show results obtained by the electric field with 
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Fig. 3. (Color online) Time dependence of hole densities dur- 
ing (t < Ti rr ) and after (t > T[ rr ) photoexcitation along the 
stripes for (a) #-RbZn (u e xt=0.4, | E ext \= 0.25) and (b) a-l 3 
(^ext=0.4, | .Eext |= 0.04). The time dependence of modulations 
in transfer integrals is also shown in (c) and (d), respectively, on 
a longer timescale. 



^ext = 0.4 and polarization parallel to the stripes. 

4. Photoinduced Melting Dynamics 

Next, we discuss photoinduced melting dynamics. Fig- 
ure 3 shows the time evolution of the hole densities 
2 — (rii) and the lattice distortions during and after the 
photoexcitation with 7V ext = 15 (Ti rr = 236 correspond- 
ing to 155 fs) which is comparable to the experimental 
pulse width. 6 The unit of time is (eV) _1 , so that t — 1520 
corresponds to 1 ps. For both salts, the COs are weakened 
by photoexcitation, as seen in Figs. 3(a) and 3(b). The 
amplitude of the electric field is chosen at | i£ ext |= 0.25 
for #-RbZn and | E ext |= 0.04 for which is near the 
critical value for melting the CO. The increment in the 
total energy per site AE is 0.098 and 0.042, respectively. 
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For #-RbZn, there are two distinct hole densities at 
t = since each stripe consists of sites with equal hole 
densities. This is also true at t > if the polarization of 
E(t) is parallel to the stripes as in Fig. 3(a), whereas the 
hole densities on the four sites in the unit cell become 
different when the polarization is perpendicular to the 
stripes. For a-I%, on the other hand, the four sites are 
already distinct at t = because of the low symmetry of 
the crystal structure. 

As for the time dependence of the hole densities, reso- 
nantly excited dynamics is clearly seen for a-ls especially 
in the oscillation of the hole densities on sites A and 
A 7 , the periods of which are nearly equal to that of the 
electric field 27r/cj e xt being 15.7. For #-RbZn, the charge 
dynamics is more complex than that for a-l3 although 
cjext is near the resonance peak of the linear absorp- 
tion spectrum. Such a difference is also visible in results 
obtained by the exact many-electron wave functions on 
small clusters of 12 sites, 43 where the origin of the dif- 
ferent behaviors is interpreted with the distribution of 
energy levels above the ground state, which is denser in 
#-RbZn than in a-ls. This is because the high symme- 
try of the undistorted structure in #-RbZn causes charge 
frustration. Therefore, a larger number of excited states 
are nearly degenerate and involved in the photoinduced 
dynamics in #-RbZn, which results in complex charge 
dynamics. 

The time dependence of the lattice distortions 
and 2/52 is shown in Figs. 3(c) and 3(d), respectively. 
Since the COs are almost completely destroyed in the 
present parameters, the displacements become zero at 
some t although they oscillate with periods longer than 
2tt/co p i 1 by phonon softening, which are much longer 
than the timescales of charge-transfer excitations. For 
^ph = 0.01, the oscillation period is about 1 ps for both 
compounds, which is comparable to the experimentally 
observed value of 0.7 ps (48 cm -1 ) in the time evolutions 
of reflectivity spectra. 6 

In order to compare the efficiency of photoinduced 
melting for #-RbZn and that for CK-I3, we show time- 
averaged hole densities as a function of AE in Figs. 4(a) 
and 4(b). Here the averages are taken after t = T- 1TT over 
the period of T irr = 27r7V ex t/^ext < t < 10irN ext /u; ext . 
There is a critical value AE C above which the horizontal 
CO disappears. The value for (9-RbZn, AE C = 0.098, is 
larger than that for CK-I3, AE C = 0.042. This is a conse- 
quence of the larger lattice stabilization energy for the 
horizontal CO in <9-RbZn than that in a-I 3 . 30-34 We no- 
tice that this inequality generally holds for any excita- 
tion frequency and for any polarization (not shown). The 
robustness of the CO in #-RbZn compared to that in a- 
I3 against the photoirradiation has been observed in the 
femtosecond spectroscopy. 6 The efficiency of the photoin- 
duced metallic state by a; e xt = 0.89 eV is evaluated as 100 
(#-RbZn) and 250 (a-13) molecules/photon. The values of 
AE C in our calculations are, however, one order of magni- 
tude larger than the experimentally estimated ones. This 
may be due to the Hartree-Fock approximation, since 
AE C becomes smaller in our recent study based on the 
exact many-electron wave functions. 43 Effects of ther- 
mal fluctuations and the experimental estimation itself, 
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Fig. 4. (Color online) Time- aver ages of hole densities after pho- 
toexcitation along the stripes for (a) #-RbZn (cj e xt=0.4) and (b) 
a-Is (d; ex t=0.4), as a function of increment in total energy per 
site AE. The time-averages of modulations in transfer integrals 
2/0 and are shown in (c) and (d), respectively. 



which generally contains an error depending on the esti- 
mation of the penetration depth, are also possible reasons 
for the discrepancy. 

For #-RbZn, the time- aver aged hole distribution be- 
comes uniform above AE C . For 0-I3, on the other hand, 
the averaged hole densities on sites A and A' become 
equal above AE C while the charge disproportionation on 
sites B and C remains for larger values of AE. This is 
reminiscent of the fact that the charge disproportiona- 
tion exists even in the high-temperature metallic phase 
of a-ls because of low-symmetry configuration of transfer 
integrals. 38 41 If we increase AE further, the hole densi- 
ties on the four sites finally merge into 0.5. 

In Figs. 4(c) and 4(d), we show the averaged lattice 
displacements which indicate that both and y^ van- 
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ish at around AE C . This is reasonable because the 
distortion is essential for stabilizing the horizontal CO 
in #-RbZn, while the distortion breaks the inversion 
symmetry that guarantees the equivalence of hole densi- 
ties on sites A and A' in a-Is. 

5. Domain Growth after Local Photoexcitation 

Recall that the mechanisms for stabilizing the COs by 
lattice distortions are different in the two salts. In 0- 
RbZn, the whole hole-rich (hole-poor) stripe is stabilized 
by strengthening (weakening) the horizontally connected 
p4 (p2) bonds. In CK-I3, the metallic phase without lattice 
distortion already possesses a hole-rich site B and a hole- 
poor site C by the relation tb2 > hi- The hole-rich site A 
and the hole-poor site A' bridged by the site B are locally 
stabilized by strengthening the bond B2' and weakening 
the bond B2. Thus, local photoexcitations would easily 
weaken the CO in CK-I3, but the CO in #-RbZn would be 
robust. 

Then, we investigate the growth of photoinduced do- 
mains with weakened CO after local photoexcitation, by 
modifying the model in eq. (2) through the introduction 
of the Peierls phase factors only on the bonds that con- 
nect four sites within a unit cell of the 12 x 12-site system. 
The phase factor is set at unity on all the other bonds. 
The time evolution is calculated by the method explained 
in §2. In Fig. 5, we show snapshots of the hole densities 
at t = 200, 400 and 600, from which those at t = 
are subtracted, (*(t)|ni|*(t)) - (*(0)|rii|*(0)), during 
and after the local photoexcitation with uj ext = 0.4 and 
iVext = 15 (Ti rr = 236). The polarization is parallel to the 
stripes. In the figures, the bonds on which the photoex- 
citation is applied are located at (7.5,7), (7.5,7.5), and 
(7.5,8). Here, the location of the bond between the sites 
at (i x ,i y ) and (j x J y ) is denoted by For 
#-RbZn (CK-I3), these bonds correspond to B-D (B-A'), 
B-A, and C-A or p3 (B2), p4 (b3), and pi (b4), respec- 
tively, in Fig. 1. In order to compare the growth dynam- 
ics, we set I £ ex t |= 0.6 for (9-RbZn and | £ ex t |= 0.32 
for CK-I3, which lead to the ratio AE/AE C ~ 0.4 for both 
compounds. Here AE C is the value obtained in §4 by the 
uniform photoexcitation, whereas AE is the increment in 
the total energy per site after the local photoexcitation. 

In Fig. 5, the growth of the domain shows anisotropy, 
which is pronounced in #-RbZn. Namely, the growth of 
the domain is faster to the direction parallel to the stripes 
than to that perpendicular to the stripes for both com- 
pounds. This is caused by the anisotropy in the transfer 
integrals which are larger for the p- and 6-bonds than for 
the c- and a-bonds in #-RbZn and a-13, respectively. For 
#-RbZn, the hole densities largely deviate from those at 
t = only on the horizontal rows near the four photoex- 
cited sites. The amount of change reaching 0.5 means 
that the CO is completely destroyed in this region. In 
the vertical direction, however, the growth is strongly 
suppressed. For a-l3, the photoinduced domain expands 
over the whole system at t = 600 although the largest 
amount of change in the hole densities is smaller than 
that in <9-RbZn. 

In order to analyze the domain growth more quantita- 
tively, we show AE/AE C as a function of | E ext | 2 in Fig. 
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Fig. 5. (Color online) Hole densities at t = 200, 400 and 600 from 
top to bottom, from which those at t = are subtracted, for 
^ext = 0.4 and iV ex t = 15. The left and right panels show 6- 
RbZn for | .Eext |= 0.6 and a-l 3 for | .Eext |= 0.32, respectively, 
where the squares indicate sites and the white lines indicate hole- 
rich stripes. The photoexcitation with polarization parallel to the 
stripes is introduced only on the bonds (7.5,7), (7.5,7.5), and 
(7.5,8), which correspond to B-D (B-A'), B-A, and C-A or p3 
(B2), p4 (b3), and pi (b4) for <9-RbZn (a-I 3 ), respectively. 
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Fig. 6. (Color online) AE/AE C as a function of | E e *t | 2 for 
#-RbZn and 01-13 after local photoexcitation with polarization 
parallel to the stripes, o; e xt = 0.4, and iV ex t = 15. 



6. For both compounds, AE/ AE C has a maximum value 
as explained in Appendix, which is below unity. Thus, 
the CO is not globally destroyed by the local photoex- 
citation with small pulse width. The maximum value is 
larger for a-l3 than for #-RbZn, which means that it is 
more difficult to weaken the CO in #-RbZn by the local 
photoexcitation with this pulse width. 

The absolute values of the differences between the hole 
densities at t = and those at t = 200, 400 and 600 are 
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Fig. 7. (Color online) The -dependence of |An^| a,t t = 
200, 400, and 600 for (a) | E e *t \= 0.3 and (b) | E e *t |= 0.8 
in the case of #-RbZn, where i y is the coordinate along the c- 
axis. The photoexcitation is along the stripes, u e xt = 0.4, and 
iVext = 15. 
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Fig. 8. (Color online) The iy-dependence of \Arii\ at t = 
200, 400, and 600 for (a) | E e *t \= 0.4 and (b) | E e *t \= 0.8 in 
the case of CK-I3, where i y is the coordinate along the a-axis. The 
photoexcitation is along the stripes, u e xt =0.4, and iV e xt = 15. 



averaged over the direction parallel to the stripes and 
denoted by \Arii\ : 

7-El<*(*)l"*l*(*)> " <*(0)|ni|*(0)>|, (12) 

i 

where L x = 12 is the number of sites along the a- and 
6-axes for #-RbZn and a-13, respectively. Thus, we can 
measure the growth of the photoinduced domain to the 
direction perpendicular to the stripes. The values |An^| 
are plotted as a function of i y in Fig. 7 for #-RbZn and 
in Fig. 8 for CK-I3, where i y is the coordinate along the 
c- and a-axes (which are perpendicular to the stripes) 
for #-RbZn and a-13, respectively, as used in Fig. 5. Al- 
though the charge order on the stripes containing the 
photoexcited bonds (7.5,7), (7.5,7.5), and (7.5,8) is weak- 
ened and the quantities |An^| at i y = 7 and 8 are in- 
creased by increasing | I2 e xt | from 0.3 to 0.8 for #-RbZn, 
the photoinduced domain remains localized near i y = 7 
and 8. In other words, the photoinduced domain hardly 
grows to the direction perpendicular to the stripes. The 
influence of the local photoexcitation is confined in a 
region near the photoexcited bonds. This property pre- 
vents AE/AE C from becoming large beyond the max- 
imum value, as shown in Fig. 6. For the domain 
growth is qualitatively different from that in #-RbZn. 
The photoinduced domain expands to the vertical direc- 
tion as shown in Fig. 8(a) for | i2 ex t |= 0.4, although the 



values of \Arti\ are relatively small. When we increase 
I i^ext |, a large region of suppressed CO (i.e., large |An^|) 
is created as shown in Fig. 8(b). It is noted that the in- 
dependence of I An^l is asymmetric for a-ls. The photoin- 
duced domain expands to the direction of decreasing i y 
more easily than to the direction of increasing i y . This 
feature is generally observed irrespective of the choice of 
photoexcited bonds because it reflects the asymmetry of 
the charge distribution that allows the ferroelectricity. 44 
Let us focus a region near the bond B2' connecting sites 
A and B that has the largest transfer integral £b2' (see 
Fig. 1). The transfer integrals t^i and £54 connecting sites 
A and C are different from ts2 an d hs connecting sites 
B and A 7 . When the local distortion at site B making 
tB2' > ^32 is weakened by photoexcitation, the hole is 
locally transferred from site A through site B to site A'. 
This corresponds to the direction that weakens the fer- 
roelectric polarization. This results in the asymmetry of 
|Arii|. It is in contrast to |An^| for #-RbZn, where the 
transfer integrals t p connecting sites A and C are the 
same as t p connecting sites B and D. 

These numerical results suggest that a macroscopic do- 
main is much more easily created in a-13 than in #-RbZn. 
They are consistent with the experimental observations 
by the time evolutions of reflectivity spectra, 6 which in- 
dicate the growth of a metallic domain only for a-ls, 
although spatial inhomogeneity is not directly observed. 
The critical slowing down observed in a-13 6 is beyond the 
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scope of the present study. If we consider it theoretically, 
we would need to treat larger system and domain sizes. 
Effects of relaxation, which are not taken into account 
here, may also be important for the dynamics on longer 
timescales. 

The different growth dynamics is caused by the dif- 
ferent mechanisms of stabilizing the COs: the CO in 6- 
RbZn is stabilized stripe by stripe through the homoge- 
neous modulation of transfer integrals along the stripes, 
while the CO in a-l3 is stabilized locally, as explained in 
§2. The growth dynamics is thus qualitatively unchanged 
even if we use different excitation frequencies and polar- 
izations in the calculations. 

6. Summary 

We have investigated the photoinduced melting dy- 
namics of COs in quasi-two-dimensional organic conduc- 
tors #-RbZn and a-ls- Although they show COs with 
similar horizontal-stripe patterns, relative importance of 
e-ph couplings and their configurations of transfer inte- 
grals in stabilizing the COs are quite different. 30 34 By 
numerically solving the time-dependent Schodinger equa- 
tion within the Hartree-Fock approximation, charge and 
lattice dynamics are obtained in the extended Peierls- 
Hubbard models during and after the oscillating electric 
field is introduced with and without spatial dependence. 
We find different photoinduced dynamics in these salts, 
which originate from different mechanisms of stabilizing 
the COs by lattice distortions. 

In the case of a spatially uniform time-dependent elec- 
tric field, we calculated time-averaged hole densities as a 
function of the increment in the total energy per site AE. 
We find a critical value AE C above which the horizontal 
CO and the lattice distortion simultaneously disappear. 
The AE C value in #-RbZn is larger than that in a-ls, 
as a consequence of the larger stabilization energy for 
the CO in #-RbZn. When the COs are destroyed, the 
hole distribution reflects the symmetry of the underlying 
crystal structure without lattice distortions. It becomes 
uniform for #-RbZn whereas the charge disproportiona- 
tion between sites B and C remains for a-ls. 

When the applied electric field is local, the time evolu- 
tion of the hole densities shows anisotropy. The influences 
of the local photoexcitation propagate more rapidly to 
the direction parallel to the stripes than to that perpen- 
dicular to the stripes. In particular, for #-RbZn a pho- 
toinduced domain hardly expands to the perpendicular 
direction. The increase of AE / AE C is limited to a small 
value when the pulse width is small, indicating that it 
is difficult to melt the CO in #-RbZn by the local pho- 
toexcitation. This is because the ^-distortion homoge- 
neously modulates the transfer integrals on the horizon- 
tally connected bonds, so that the charge correlation in 
each stripe (between stripes) is strong (weak). For a-ls, 
on the other hand, the photoinduced domain expands 
in the plane. This comes from the fact that each hole- 
rich bond B2' between sites A and B in a-ls is locally 
stabilized by the ^-distortion. 

The obtained results are qualitatively consistent with 
the experimental observations, which indicate that the 
CO in #-RbZn is more stable against photoexcitation 



than in a-ls: a macroscopic metallic domain is generated 
in 0-I3 whereas the CO only locally melts in #-RbZn. 
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Appendix: Maximum of Increment in Total En- 
ergy 

Figure 6 shows that AE/AE C begins to decrease at 
around | £ ext |= 0.7 (| E ext | 2 = 0.5). This behavior 
comes as a consequence of the periodicity of the en- 
ergy band in momentum space. The electric field shifts 
the momenta of electrons. For simplicity, let us consider 
the noninteracting electrons in one-dimension, where the 
shift is estimated by the equation of motion, 

dk 

— = E ext sincj ext t. (A-l) 
at 

Here k is the wave number and we set e = a = h = 1 with 
a being the lattice constant. By integrating the above 
equation from t = to t = 7r/uj ext that is one half of 
the period of the oscillating electric field, we obtain the 
momentum shift \Ak\ as, 



^ext 

As E ext increases, \Ak\ and AE increase. When \Ak\ ~ 
7r, however, some electrons accelerated by the electric 
field cross the boundary of the Brillouin zone so that \Ak\ 
is no longer proportional to E ext . If we increase E ext fur- 
ther, it is expected that the total energy is not efficiently 
increased any more and even decreases because more 
electrons cross the zone boundary. This results in the 
decrease of AE / AE C as shown in Fig. 6. For c^ext = 0.4, 
| A A: | ~ 7r corresponds to E ext ~ 0.6, which is consistent 
with the above value, | E ext |= 0.7. 
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